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Maximally entangling tripartite protocols for Josephson phase qubits 
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We introduce a suit of simple entangling protocols for generating tripartite GHZ and W states 
in systems with anisotropic exchange interaction g {XX + YY) + gZZ. An interesting example 
is provided by macroscopic entanglement in Josephson phase qubits with capacitive {g = 0) and 
inductive (0 < \g/g\ < 0.1) couplings. 
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I. INTRODUCTION 



Superconducting circuits with Josephson junctions have attracted considerable attention as promising candidates 
for scalable solid-state quantum computing architectures. The story began in the early 1980's, when Tony Leggett 
made a remarkable prediction that under certain experimental conditions the macroscopic variables describing such 
circuits could exhibit a characteristically quantum behavior [1| . Several years later such behavior was unambiguously 
observed in a series of tunneling experiments by Devoret et al. Q], Martinis et al. 3, and Clarke et al. [4|. It was 
r^ I eventually realized that due to their intrinsic anharmonicity, the ease of manipulation, and relatively long coherence 
Cd ' times Q, the metastable macroscopic quantum states of the junctions could be used as the states of the qubits. That 
? , . idea had recently been supported by successful experimental demonstrations of Rabi oscillations 6], high-fidelity state 
^^- preparation and measurement [3, @, M, UQ, \lM, \lM ' ^'^'^ various logic gate operations [8, 9., ^0, JJ,, 1S\ . Further progress 
] in developing a workable quantum computer will depend on the architecture's ability to generate various multiqubit 
—i • entangled states that form the basis for many important information processing algorithms il4|. 
^ I In this paper we develop several single-step entangling protocols suitable for generating maximally entangled quan- 

^^ . turn states in tripartite systems with pair-wise coupling g {XX + YY) + gZZ. We base our approach on the idea 
that implementing symmetric states may conveniently be done by symmetrical control of all the qubits in the system. 
This bears a resemblance to approaches routinely used in digital electronics: while an arbitrary gate (for example, a 
3-bit gate) can be made from a collection of NAND gates, it is often convenient to use more complicated designs with 
three input logic gates to make the needed gate faster and/or smaller. 

The protocols developed in this paper may be directly applied to virtually any of the currently known superconduct- 
^^ ' ing qubit architectures, two of which will be mentioned here. The first architecture is based on capacitively coupled 
current-biased (CBJJ) Josephson junctions [1^, [iSJ, |l5[ whose dynamics is governed by the circuit Hamiltonian 



X ■ Hi = {pI+pI + KP1P2) /2m + (?i/2e) V [-/q cos </>^ - /,0,] , (1) 



^2 



with Pi — m(f)i, m = {h/2e) (0 + Cint), k = 2Cint/(C -I- Cint). The other architecture involves inductively coupled 
flux-biased (FBJJ) junctions [la]. It is described by the Hamiltonian (for small T, see Ref. [131 for details) 

2 
H2 = {pI +pI) /2m + {h/2e) ^ [-/q coscj), + {eEo/h) (0, - 27r$,/$3c)'] + TEq (0i - 2^$,/$3c) (02 - 27r$,/$3c) , 



i=l 



(2) 
with p, = m0j, m = {h/2efC, wq = I/VLC, Eq = h^io^/2Ec, Ec = (2e)V2C, $sc = h/2e, T == M/L. When 
reduced to computational subspace, in the rotating wave approximation, these Hamiltonians become 

Hkwa = (1/2) [fli ■di+n2-a2+g {alal + a>^) + gala^] , (3) 
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with ^ = (momentum- momentum coupling) and < g/g < 0.1 [T^ (position-position coupling), respectively. 
For typical superconducting qubits, the RWA requirements are well satisfied: the level splittings are usually around 
uj/2'n: ~ 10 GHz and the Rabi frequencies fi needed to implement various logic gates [13, lly] are on the order of the 
coupling constant g/ln < 100 MHz. Thus, il/uj ~ 10~^, as required. 



II. THE GHZ PROTOCOL 



A. Triangular coupling scheme 



In the rotating frame in the absence of coupling, the computational basis states |000), |001), |010), |100), |011), 
|101), |110), |111) have the same effective energy Ecff = 0. The pair-wise couphng, 
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35/2/ 



(empty matrix elements are zero), partially lifts the degeneracy, which results in the energy spectrum 

£;int = {35/2, 35/2, 2g - 5/2, 2g - ~g/2, - (<? -f ~g/2) , - (5 -f ~g/2) , - {g + ~g/2) , - (5 -f ~g/2)} , (5) 

with the corresponding H-eigenbasis 

HgHZ 0Ww0Hrest = {|000) |111)}0{|W) |W')}0 {|*l) ® |*l') © |*2) © |*2')} , (6) 

where 

|W) = (|100) + |010) + |001))/\/3, |W') = (|011) + |101) + |110))/V3, 

\^,) - (|100> - |010))/V2, \^,') = (1011) - |101»/V2, 

1*2) = (|100) + |010)-2|001))/\/6, |*2') = (|011) + |101)-2|110))/V6. (7) 

Since the coupling does not cause transitions within each of the degenerate subspaces (nor does it cause transitions 
between different such subspaces), it is impossible to generate the |GHZ) = (|000) + |lll))/-\/2 state [20] from the 
ground state |000) by direct application of i/int- Instead, we must first bring the |000) state out of the TYghz subspace 
by, for example, subjecting it to a local rotation i?i in such a way as to produce a state ji/') that has both |000) 
and I 111) components. That is only possible if all one-qubit amplitudes ai,...,(3z in the resulting product state 
\'4>) = i?i|000) = (ai|0) + (3x\l)) (a2|0) -|-/32|1)) (aajO) -|-/?3|1)) are chosen to be nonzero, which means that in the 
computational basis the state |^) will have eight nonzero components. 

We now notice that in the 7i-basis, the three-qubit rotations 
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-VSsc^ -s(l-3c2) c(l-3s2)/ 



where Yg ' = exp (— i6'o'*;'/2), Xg ' = exp (— i0cr^'/2), k — 1,2,3, are block-diagonal, with c = cos(0/2) and s 
sin(6'/2). For 6 — tt/2, the corresponding 4x4 blocks acting on the TYghz © Ti-w subspace are 

/ 1 i -iVS -VS\ / 1 -1 



3^(4x4) _ ^ 
^7r/2 - -^ 



1 



i — i\/3 — \/3\ 

i 1 — v3 —iVS 

-iVS -\/3 -1 -i 

V-\/3 -iV3 -i -1 / 



V(4x4) _ ' 
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V8 



-V3 Vi\ 
1 1 V3 Vs 
V3 \/3 -1-1 
VV3 -Vi 1 -1/ 



(9) 



This shows that y^/2 provides a convenient choice for i?i . We may thus start by generating the so-called symmetric 
state, 

IV').y„, = i^V2|000> = (1/2) (|GHZ> + v/372 (|W> + |W'))) e Hghz 0Hw. (10) 

The entanglement is then performed by acting on |'0)svin "with C^int = exp (— ii/jnt^), thus inducing a phase difference 
between the GHZ and W+W components (this step works only for g 7^ g, see Section ITV)) . 

/2|000) = (e-^72)(|GHZ)+e-^^V372 (|W> + |W'))), a^{3g/2)t, S = 2{g-g)t. (11) 



CAnti^TT, 

To transform to the desired GHZ state, we first diagonalize the X^,^ ' and Y^,^ ' operators to get the unimodular 
spectra 



Xx = (-e'(''/'*),-e"'(''/'*\e"'(''/'*\e*(''/'*^| , Xy = |-e"'(''/'*\ -e'^''/''', e*(''/"*\ e^'^'^/^^l 



(12) 



corresponding to the X- and J-eigenbases, X = (|Xi) jXa) IX3) IX4)) = y}^^^\ y = {\Yi) lY^) iFs) 1^)) = 



X^ /I ' , which are formed by the columns of y ,^ and X^ /I ' . Using the A'-basis, we notice that both states 



IGHZ) = ^il±^^Mil , c/,,i;/,|ooo) 



'-7r/2 



l + 3e 



-4(5 



\Xi) 



1-e 



-i5 



-V3IX4 



(13) 



belong to the same two-dimensional (nondegenerate) A'-subspace spanned by \Xi) © 1^4) ■ Therefore, by performing 
an additional X^/2 rotation we can transform C/mty7r/2|000) to 



^V2t^int>;/2|000) = e-^"e'(-/4)|(jjjZ), 



(14) 



provided the entangling time is set to give \6\ = tt, or, icHZ = 7r/2|(7 — g\. Any other GHZ state (|000) -|- e"^|lll))/V2 
can be made out of the "standard" GHZ state by a Z-rotation applied to one of the qubits, as usual. 

The protocol may be compared to controUed-NOT logic gate implementations [13, Il3l that used various sequences 
CNOT = e-*(''/'*)i?2t/cNOT^i, with (entanghng) times icNOT = TTT/2g, 1<T < 1.6. Thus, for g = 0, the entanghng 
operation proposed here will be of same duration as the fastest possible CNOT. 

We conclude this section by noting that in its present form the GHZ protocol cannot be used to generate the W 
state. This can be seen by writing |W) = (\/3(|Xi) -I- IX2)) - (IX3) + IX4))) /V&, which shows that our XUintY 
sequence does not result in a W since the final X^/2 rotation cannot eliminate the IX2) and {X^} components. Also, 



|W> = (V3 (*|yi) - 112)) - (1^3) - ^\Y,))) /VS, 



(15) 



and 



I ^ — '^p-''' \/3 (1 + e 

Y./2UintY,/2\000) = e-^" ^ m,) - \Y2)) - ^ ^ 



-i5\ 



(l>3>-*in» /V8, 



(16) 



and thus no choice of S will work for the YUintY sequence either. 



B. Linear coupling scheme 



In the case of linear coupling, say, 1 <-> 2 and 2 <-> 3, the interaction Hamiltonian is given by 
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{5, 3, e(+\ e^+\ e^-\ e^-\ 0, o} , e^ - ±^2g^ + {~g/2f - ~g/2, 



(17) 



with eigenbasis 

1000), |111), 

|W)^+^ = C(+) (|001) + (e(+V.9)|010) + |001)) , jW^^^ = C(+) (|011) + (e(+V,g)|101) + |110)) 

|W)'"^ = d-) (|001) + (e(-V5)|010) + |001)) , |W')*"' = C^'^ (|011) + (e(-V5)|101) + |110)) 
\^) = (1001) - |100))/n/2, I*') = (1011) - |110))/n/2, 
where C^ ' are normaUzing constants. We have, 



|W) =aW|w)(+^+a(-)|w)'"\ aw = 



and similarly for |W'). Our GHZ sequence then leads to the entangled state 



A(-) = 



C/int>^V2|000) = (e-'"/2) (|GHZ) + ^3/2 (^e-'^'^' A'^+I [|W)(+^ + |W') 
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(18) 
(19) 



(20) 

with a ~ gt, (5'^*' = (e'^^'— g)i. Since i > 0, in order for the X^/2 post-rotation to give a GHZ, we must restrict coupling 
to g = and set the entangling time to <ghz ~ ^/v2|<?|- An alternative GHZ implementation for superconducting 
qubit systems with capacitive coupling has recently been considered in Refs. [2ll. |22|. There, individual qubits were 
conditionally operated upon one at a time. 



III. THE W PROTOCOL 

We now turn to the W protocol. Eq. ([TB)) suggests that control sequence YUinty may still give a W, provided a 
proper adjustment of i\Yi) — {¥2} and jYa) — i|Ki) amplitudes is made by a physically acceptable change of system's 
Hamiltonian. In the context of Josephson phase qubits such modification can be achieved by adding local Rabi term(s) 
to Hint, for instance, 

/3g n n n \ 
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n 2.9 ~g 2g n fl 

n 2q 2g -g n fl 
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Hl^,, = {n/2){al+al+al)+Hi, 



(21) 



The energy spectrum then becomes E^^ 



g/2 ± n/2, x^"^^ = ^/{g~ g)^± {g -g)n + n'^. The (first two) eigenvectors are 



l*r,^> - ciV 



-1- 



|e(+) ± ^(+), e(-) ± ^i-)^ _,(+), _£(+), _£(-), _e(-)}, ^ith e(±) =g + 

(22) 
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^(+) i._ 



with normalizing constants C^ , fc = 1, 2. After some algebra we find: 



f/Sti;/2|000) = e-'"/{4V2x'^+'>) [{A/n) {te'^-^/^^\Yi) - e-^^'/^^lFa)) + {V3B/n) (e-^f'/^^lFs) - Je'^^/^^l 



Y. 



where 



A = (^g-g + n + x(+)) [g-~g + 2n- x^+^) - e"'* (^g - ~g + n - x^+^) (g - g + 20 + x(+^) 
B = (^g-g + n + x^) (.9 - ~g - x'+') - e"*' (.9 - g + - x(+)) (5 - .9 + X<+^) , 



(23) 



(24) 



and 



(' 



:( + ) 



C~^') t, S = —2x''h- It is straightforward to verify that additional Yt^/2 rotation applied to this 



state produces a W (see Eqs. P^. (ITS|) '). 

rv2C/St>^V2|000) = hsgn (5 - g)] e~'"\W), 
provided we set i-yv — T^/V^lg — g\, ^ = —{g — ff)/2- 



(25) 



IV. ADDENDUM: ISOTROPIC HEISENBERG EXCHANGE g{XX + YY + ZZ) 

Maximally entangling protocols introduced in previous sections are singular in the limit g —^ g, which corresponds 
to the isotropic Heisenberg exchange interaction. Even though this limit is not met in superconducting qubits, for 
completeness, we briefly discuss it here. 

It is obvious that when g = g, the symmetric state y^/2|000) is an eigenstate of the interaction Hamiltonian. 
Consequently, the Heisenberg exchange does not cause transitions out of it, making the gate time divergent. To 
perform single-step entanglement we break the symmetry of local rotations. For example, the GHZ state can be 

generated by e-'"|GHZ) = e-*('^/2)'-'e~*(^/3H'^i-'^«)[/i„te-''(^/i^n5'^i+'^2-3'^«)e-»W2)'^'|000), with a = -n/2, toaz = 
(2/3) X (Tr/2g). To generate the W state, we generalize Neeley's fast implementation for triangular g(XX + YY) 
coupling [23] (c/. [2J]) to arbitrary coupling g{XX + YY)+gZZ, including the Heisenberg exchange g = g: e~*"|W) = 
g+»(V3)<T^y.^^g-»(V2)<Tj|QQQ^^ with a = (55 - 2g)7r/185, iw = (4/9) x (7r/25). 

V. CONCLUSION 

In summary, we have developed several single-step symmetric implementations for generating maximally entangled 
tripartite quantum states in systems with anisotropic exchange interaction, which are directly applicable to supercon- 
ducting qubit architectures. In the GHZ case, both triangular and linear coupling schemes have been analyzed. In 
the isotropic limit, our implementations exhibit singularities that can be removed by breaking the symmetry of the 
local pulses. 
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